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Abstract 

Wc define the Euler number of a bipartite graph on n vertices to be the number of labelings of 
the vertices with 1, . . . , n such that the vertices alternate in being local maxima and local minima. 
We reformulate the problem of computing the Euler number of certain subgraphs of the Cartesian 
product of a graph G with the path P m in terms of self adjoint operators. The asymptotic expansion 
of the Euler number is given in terms of the eigenvalues of the associated operator. For two classes 
of graphs, the comb graphs and the Cartesian product P2^P m , we numerically solve the eigenvalue 
problem. 



1 Introduction 



Let G = (V,E) be a bipartite graph on n vertices with the vertex decomposition V = V\ U V2, that 
is, each edge in G has one vertex in V\ and the other in V2. An alternating labeling n is a bijection 
7r : V — > {1, . . . ,n} such that for two adjacent vertices u £ V\ and v E V2 we have that ir(u) < tt(v). 
Another way to phrase this condition is that every vertex u in V\ is a local minimum of the bijection ir 
and every vertex v in V2 is a local maximum. Define the Euler number £{G) to be the number of 
alternating labelings tt of the vertices of the graph. 

Two examples are in order. First, for the path P n on n vertices the Euler number £{P n ) is the 
number of alternating permutations, that is, the classical Euler number E n . Second, for a cycle C n of 
even length n we have that £(C n ) = n/2 ■ E n -\. Since there are n/2 possible positions for the largest 
label n, the labeling reduces to the path P n -\. 

Observe that we cannot drop the condition that the graph G is bipartite, since the labeling can 
not alternate along an odd cycle. Alternatively, for a non-bipartite graph G let £ (G) = 0. Observe 
that the definition of the Euler number is independent of the order of V\ and V2. We also have the 
trivial lower bound 

\V 1 \\-\V 2 \\<£(G), (1.1) 

by assigning V2 the | V2 1 largest labels. Equality in (jl.ip is only obtained for the complete bipartite 
graphs. Moreover, extending the classic "Multiplication Theorem" due to MacMahon [4, Article 159], 
for the disjoint union of two graphs G and H we have 

£(G U H) = ( m + ^ ■ £(G) ■ £(H), (1.2) 
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Figure 1: The comb graph Comb m as the product P 2 0^ 1 yP 1 



where G and H have m respectively n vertices. 

Our interest is to study subgraphs of the Cartesian product of two graphs. For graphs G and H, 
and S a subset of vertices of G, define the product Go s H as the graph on the vertex set V{G) x V(H) 
where two vertices (u, u') and (v, v') are adjacent in GO^H if either u = v G S and u' is adjacent to 
v', or u is adjacent to v and v! = v'. The Cartesian product GOH is obtained as a special case of the 
product GU S H with S = V{G). 

For the general problem we obtain the following asymptotics. 

Theorem 1.1. Let G be a bipartite graph on n vertices and S a non-empty subset of the vertices of 
the graph G. Then there exist three positive real numbers X, fx and c such that A > fi and 

£(GO s P m ) , m _i . m m -\\ 

— — = c-A +U(fi ) asm — > oo. 

(m ■ n)\ 



2 The self adjoint operator 

For a bipartite graph G = (V, E) on n vertices define two subsets X and Y of the n-dimensional unit 
cube [0, l] v in n-dimensional Euclidean space R v by 

X = {xe[0, 1} V : x u + x v < 1 for {u,v} G E}, 

Y = {x G [0, l] V : x u < x v for u G V\, v G V 2 , {u, v} G E}. 

Lemma 2.1. The two subsets X and Y have the same volume, which is given by the Euler number 
of the graph G divided by n\. 



Proof. By reflecting the set Y over all of the hyperplanes of the form x v = 1/2 where v G V 2 we obtain 
the set X. Hence their volumes agree. By cutting the n-dimensional cube with the hyperplanes 
x u = x v for all u, v G V we obtain n! simplices of the same volume. Each simplex corresponds to a 
permutation by reading the order of the coordinates of a point in its interior. The set Y is the union 
of a subcollection of these simplices corresponding to an alternating labeling of the graph G. □ 



Let the 0, 1-function x be defined on the set X x X by 

x(x, y) = 



1 if x v + y v < 1 for all d£5, 
otherwise. 



Definition 2.2. Define the operator T on L 2 (X) by 

T[f](x)= [ x(x,y)-f(y)dy. 
Jyex 

Since % ls symmetric, that is, x(x,y) = x{v-,x), the operator T is a self-adjoint Hilbert-Schmidt 
operator. Thus we conclude that the spectrum of T is real and discrete with as the only possible 
accumulation point. Furthermore, all the eigenvalues and eigenfunctions of the operator are real. 
Since < x(x,y) < 1, the eigenvalues A lie in the closed interval [—1,1]. Hence, there is a largest 
eigenvalue in absolute value. Moreover, the eigenfunctions form a complete orthogonal set. 

Let 1 denote the constant function with value 1 on set Y. Now we have 
Proposition 2.3. For a bipartite graph G on n vertices and S a subset of the vertices of the graph G, 

£{Gu s P m 



(m • n)\ 



<l,T m - 1 [l]>. 



Proof. Expanding the inner product and each of the m—1 applications of the operator T, we have 
that 



J xi 

J XI Jx-2 



X\ JX 2 



X{xi,x 2 ) ■ x(x2,x 3 ) ■ ■ ■ x(x m -i,x m ) dx m ■ ■ ■ dx 2 dxi. (2.1) 



Let x*i be the vector (x V i) V £y Then the integral in equation (12. 1|) is over all of the m ■ n variables x v j 
with the boundary condition that (i) < x V; i < 1, (ii) x v> i + x v ^ + \ <lforl<i<m— 1 and v £ S, 
and (hi) x U) i + x v ^ < 1 for {u, v} an edge in G. These inequalities describe exactly the set Xcn s p m 
and hence the integral is given by the ratio £(Go s P m )/(m ■ n)\. □ 

Theorem 2.4. Let G be a bipartite graph on n vertices and S a non-empty subset of the vertices of 
the graph G. Then we have 

£{GU s P m ) (y fc ,i) 2 ^ m _ 1 

(m-n)l ^ ||(Pfc|| 2 k ' 

where the eigenvalues of the operator T are {\k}k>o and fk *s the eigenfunction associated to the 
eigenvalue A&. 



Proof. Expand the function 1 in terms of eigenfunctions: 

\- (<p k ,l) 

\\Vk\\ 
fc>i ll ^' c|1 

Apply T m_1 and take the inner product with 1 and the result follows. □ 
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When the set S is empty, Theorem 12.41 is trivial. In that case Go l j l P m is the disjoint union of m 
copies of G. Using (|1.2p we have that 

£(Ga 9 P m ) = f£(G)\ m 
(m ■ n)l \ n! J 

Example 2.5. When the graph G consists of a singleton vertex and S consists of this vertex, then 
the product Go s P m is exactly the path on m vertices, and its Euler number is the classical E m . In 
this case the operator T is given by 

T[f](x)= f X f(z)dz. 
Jo 

This operator has eigenvalues Afc = 2/(ir • k) where k = ...,—7,-3,1,5,9,... and eigenfunctions 
ifk = cos(x/Afc). Calculating = Xk and ||</?fc|| 2 = = 1/2 we obtain the following 

classical asymptotic expansion for the Euler number 

(2 \ m+l 
—k) 

= 2- ml- Y, . f AS\ 

j odd 

where j = \k\, that is, k = (-l)O'- 1 )/ 2 • j. See [21 Section 4]. 



(m+l) 



Let W be the subspace of L 2 (X) consisting of the functions only depending on the variables x u 
where u belongs to S. In the case when S is the vertex set of the graph G the space W is L 2 (X). The 
following result applies to the case when S is strictly contained in the vertex set of G. 

Proposition 2.6. The image of the operator T is contained in subspace W. Hence all the eigenfunc- 
tions associated to non-zero eigenvalues belong to W . 



Proof. For a vertex v not in the set S, observe that the function x(af, y) does not depend on the 
variable x v . Hence when integrating x(x,y) ■ f(y) over all y £ X the resulting function T[f] does not 
depend on x v , that is, T[f] belongs to the space W. The second statement follows from the defining 
relation for eigenfunctions. □ 



The Frobenius-Perron result applies to matrices, that is, linear operators on a finite-dimensional 
vector space. An operator version of Frobenius-Perron was discovered by Krein and Rutman [3j. We 
present a specialized version of their result. Let Z be a measurable space. Recall that two functions 
in L 2 (Z) are considered the same if they differ on a set of measure 0. We call a function / € L 2 (Z) 
non-negative if f(x) > for almost all x £ Z. Similarly, we call the function / positive if f{x) > for 
almost all x £ Z. An operator M on L 2 (Z) is positivity improving if for all non-negative but non-zero 
functions / the function M[f] is positive. 

Theorem 2.7 (Krem-Rutman). Let M be an operator L 2 (Z) such that there is a positive integer k 
so that M k is positivity improving. Then the largest eigenvalue X (in modulus ) of M is real, positive 
and simple. Moreover, the associated eigenfunction if is a positive function on Z. 
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m 


E 


£ (Comb TO ) 


~ £(Comb m ) 


£(P 2 aP m ) 


~£(P 2 OP m ) 


1 


1 


1 


0.99379166 


1 


0.98451741 


2 


1 


5 


4.99961911 


4 


4.002193 


3 


2 


66 


65.9990972 


44 


43.99713 


4 


5 


1613 


1612.99965 


896 


896.0018 


5 


16 


63480 


63480.0072 


29392 


29391.93 


6 


61 


3662697 


3.66269757 x 10 6 


1413792 


1.413789 x 10 6 


7 


272 


291407424 


2.91407470 x 10 s 


93770800 


9.377064 x 10 7 


8 


1385 


30572578425 


3.05725832 x 10 10 


8201380224 


8.201366 x 10 9 


9 


7936 


4089549416832 


4.08955006 x 10 12 


914570667792 


9.145691 x 10 11 


10 


50521 


679329771871725 


6.79329879 x 10 14 


126651310675680 


1.266511 x 10 14 



Table 1: Table of the Euler numbers, the number of alternating 2 x m arrays, comb graph Comb m , 
and their numerical approximations, denoted by ~ £(Comb m ) and ~ £ {P^Pm). 

Applying Krem-Rutman to our operator T, we have 

Proposition 2.8. The operator T 2 is positivity improving. The largest eigenvalue (in absolute value) 
A of the operator T is real, positive and simple. Furthermore, the associated eigenfunction (p is positive. 

Proof. Let / be a non-negative, non-zero function in L 2 (X). By the definition of the operator T we 
have that in a neighborhood of the function T[f] has a positive support. By applying the operator T 
again we obtain that every point in the interior of Y takes a positive value in the function T 2 [/]. The 
remainder of the proposition follows from Krem-Rutman. □ 

Proof of Theorem li.il By letting A be the largest eigenvalue of the operator T and letting [i be a 
bound on the next largest eigenvalue such that A > fx, the result follows. □ 



3 The comb graph 

We now turn our attention to the comb graph. See Figure [TJ Recall that the comb graph is defined 
by the product P2^{i}P m - In this case the space X is the triangle 

X = {(x,y) : x,y>0,x + y<l}. 

However, following Proposition 12 . 61 in order to find the eigenvalue and eigenfunctions of T it is enough 
to consider the subspace W of L 2 (X) consisting of functions depending only on the variable x. Observe 
that W inherits the inner product 

{f,g)w = (! ~ x ) ■ /( x ) • 9{x) dx. 
Jo 

Moreover, the operator T is given by 

T[f](x)= f 1 X (l-z).f(z)dz. 
Jo 
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x 10" 



Figure 2: The difference /(1/2) — g(l/2) found by solving the system of 
ODE in (|3.3p with a given value of A. The roots of this plot correspond to 
eigenvalues. The lower plot is a magnification of the center domain. 



The next step is to find all of the eigenvalues and eigenfunctions of the operator T, that is, functions 
f(x) so that 



l-x 



X-f(x) = T[f)(x) = / (l-z)-f(z)dz. 
Jo 



(3.1) 



We convert this integral equation into a differential equation by differentiating to obtain 

A • f'(x) = -x ■ /(l - x). (3.2) 
To convert this into an ordinary differential equation (ODE), we define g{x) = f(l — x) and thus 



-x/X 
{l-x)/X 



Together with the boundary conditions 

/(0) = 1, 5 (0) = 0, /(l/2)= 0(1/2), 



(3.3) 



(3.4) 



which come from the integral equation (|3.ip and the algebraic relationship between / and g, this linear 
system is equivalent to the original integral equation. The condition that /(0) = 1 is our choice of 
normalization for the eigenfunctions. The only solution which has /(0) = g(0) = is identically zero, 
thus this normalization is valid. We proceed to solve this ODE numerically. 
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1=0.437141... ?i=-0.0943304... 




Figure 3: The eigenfunctions associated with the four largest (in absolute 
value) eigenvalues. The function f(x), in [0,1] is found by using f(x) in 
[1/2,1] and g(l - x) in [0,1/2]. 



First, we solve the system from to 1/2 for various values of A and find the difference /(1/2) — 
g(l/2). This is plotted in Figure [2j The roots of this plot correspond to eigenvalues and we find them 
numerically. For each eigenvalue A the functions / and g are found and f(x) over the unit interval 
is reconstructed. The eigenfunctions corresponding to the largest (in absolute value) eigenvalues are 
plotted in Figure El The resulting norms and. constants c n are tabulated in Tabled 

A </(*),!) ||/(x)|| 2 



0.437141117 0.437141151 0.398916677 0.479028320 

-0.094330445 0.094331326 0.690741849 0.012882380 

0.053662538 0.053688775 0.829794009 0.003473735 

-0.037528586 -0.037546864 0.932757330 0.001511397 

Table 2: The values of A, (f(x),l) w , \\f(x)\\yy, and c n for the first 
four eigenfunctions shown in Figure [3j The constant c n is the ratio 
(f( x ), !)w for tne n * h eigenfunction. 



The Euler numbers for the Comb n graphs are calculated from the numerical approximation for 
A n and c n using the first four terms in the series in Theorem 12. 4[ They are tabulated in the fourth 
column of Table [TJ 
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4 Alternating 2 by m arrays 



The Euler number of the Cartesian product of two paths P m and P n counts the number of alternating 
m by n arrays. That is, the number of assignments of the integers 1, 2, . . . , m ■ n to an m by n array 
such that each entry is a local maximum or a local minimum. Hence, if i + j is even then the entry 
Ojj should be less than the four adjacent entries ai-ij, flj+ij, ciij-i, Similarly, if i + j is odd 

then the entry djj should be larger than the four adjacent entries. 

In the following we study the number of alternating 2 by m arrays, that is, the Euler number of 
the graph P^^Pm- The graph G is the path on two vertices P2 and 5 = {1, 2}. As before, the space 
X is the triangle 

X = {(x,y) : x,y>0,x + y<l}. 
Observe that the operator T has the form 

T[f](x,y)= I f(z,w)dzdw, 
JR 

where R is the region described by the inequalities < z < 1 — x, < w < 1 — y and z + w < 1. Since 
x + y < 1, equivalently (1 — x) + (1 — y) > 1, the inequality z + w < 1 cuts off a triangle from the 
rectangle [0, 1 — x\ x [0, 1 — y]. Hence we have 

T[f](x,y) = f V f 1 V f(z,w)dwdz+ f X f 1 " f(z,w)dwdz (4.1) 

JO JO Jy Jo 

rx ri—x r^-—y rl—W 

= f{z,w)dzdw+ j I f(z,w)dzdw. (4.2) 

J0 JO Jx JO 

In order to study this operator T, it will be easier to work in a different space. Let U be the space 
of functions g(x) on the interval [0, 1] that satisfy the inequality 

/ (g(x) + g(y)) ■ (g(x) + g(y)) dxdy < 00. 
Jx 

Enrich the space U with the following inner product 

(9,h) u = / {g{x)+g{y))-{h{x) + h{y))dxdy. 
Jx 

Define L to be the linear map L : U — ► L 2 (X) defined by L[g](x,y) = g{x) + g(y). Observe that 
the map L preserves the inner product, that is, L is an isometry. Moreover, L is an injective map. 
Furthermore, define the operator T on U by 

f*X P 1 — -X 

T[g]{x) = (1 - x) ■ / g{s)ds+ / (1 - s) ■ g(s) ds. (4.3) 

Jo Jx 

The reason why we denote this operator also by T will be clear from the next proposition. 
Proposition 4.1. The isometry L and the operator T commute, that is, T o L = L o T. 
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Proof. Apply the original operator T to the function L[g](x,y) = g(x) + g{y). We do this by applying 
equation (|4.ip to g(x) and applying equation f|4.2j) to g(y). We then have 

nl—y rl-x rl-Z 

g{z) dwdz+ I I g(z) dw dz 
Jy JO 

g(w)dzdw + / / g(w)dzdw 
Jx Jo 
ry ri—x 
= (1 - y) • / g(z) dz+ (1 - z) ■ g(z) dz 

JO Jy 

nx rl-y 

+ (1 — x) ■ j g{w) dw + I (1 — w) ■ g(w) dw 

JO ' Jx 

= L[T[g]](x,y), 

where we used the fact that f + = ' + /?• □ 



Hence we have that 

(1,^[1]) L2{X) = \. {1,^11])^ 

since L[l/2 • 1] = 1. Thus it is enough to solve the eigenvalue problem A • g(x) = T[g](x) in U for 
non-zero A, that is, 

A • g{x) = (1 - x) ■ r g(s) ds + / *(1 - s) ■ g(s) ds. (4.4) 

Jo Jx 

Differentiate to obtain 

I-X 

A • g'(x) = — g(s) ds — x ■ g{\ — x). (4-5) 
Jo 

Differentiate again 

A • g"{x) = -g{x) - g(l - x) + x ■ g'(l - x). (4.6) 
As before, we convert this into a linear system of differential equations by defining h(x) = g(l — x). 
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y 


/ 








1 







( 9 \ 


h 















1 




h 


g' 






-1/A 


-1/A 





-x/X 




9' 






V 


-1/A 


-1/A 


(l-x)/A 





/ 


\h> J 



(4.7) 



We solve this system differential system numerically on the interval [1/2, 1]. To find boundary condi- 
tions we set x = 1/2 in equations (j4.4|) and (14. 5j) . 



r i/2 

A - 5 (l/2) = 1/2- / g(s)ds, (4.8) 
Jo 

f l/2 

\-g'(l/2) = - g(s)ds- 1/2 -g(l/2). (4.9) 
Jo 

Observe that 5(1/2) = h(l/2) = 0, gSD and imply that g'{l/2) = h'(l/2) = and therefore 
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o 




2 






1 

















-1 






-2 






-3 





-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 




Figure 4: The value of h'(l) found by solving the system of ODE's in (|4.7p 
with a given value of A. The roots of this plot correspond to eigenvalues. 
The lower plot is a magnification of the center domain. 

corresponds to the zero solution of (|4.7H . Since we are looking for the non-zero solution, we normalize 
such that 

5 (l/2) = Ml/2) = A/2 + 0, (4.10) 

This gives us two conditions at x = 1/2, and also implies that J^ 2 g(x) dx = A 2 . Combined with (|4.9p 
this gives two more conditions 

g'{\/2) = -h'{\/2) = -A - \. (4.11) 

Thus given the parameter A we can solve the system from s = 1/2 to x = 1. At x = 1 however, the 
integral equation (|4.5|) yields another constraint: 

h'(l) = 0. (4.12) 

Equivalently, looking at (|4.4p (and remembering that g{x) = h(l — x)) for x = 0, 1 we find that 

h(l) = g(0) = -g(l). (4.13) 

This is only satisfied for a discrete set of eigenvalues A. To find this set we (numerically) solve the 
ODE starting from x = 1/2 using the initial conditions given by (|4.10[ I4.11|) and search values of A 
that give h'(l) = 0. Figure H] shows the value of h'(l) for various values of A. 
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A numerical root-finding algorithm finds the first few roots, that is, eigenvalues A. The associated 
eigenfunctions are shown in FigureO Finally, to find c n we must evaluate {g, 1} U and (g,g) u . This is 
again done numerically. The results are shown in Table [3j 

^=0.364426... ^=0.06401 91... 




^=-0.061 41 98... A,=0.0327004... 

1 i ' 1 . ' 1 2 i ' ' 1 



0.5 



-0.5 



0.2 



0.4 



0.6 



0.8 




Figure 5: The eigenfunctions associated with the four largest (in absolute 
value) eigenvalues. The whole function g(x), from to 1 is found by using 
g(x) between 1/2 and 1 and h(l — x) between and 1/2. 



(<?(*),!> 



A 

0.364425573038 
0.064019105418 
-0.06141983509 
0.03270035262 



0.59039705924381 
-0.05366486422899 
0.04799387821016 
0.02267022514715 



l|g(x)|j 2 

0.75905252379149 
0.29041608589489 
0.10956972701418 
0.24295945072711 



0.45921550437989 
0.00991652250888 
0.02102234265267 
0.00211532873771 



Table 3: The values of A, (g(x), 1}^, ||g(x)||jy, and c n for the first four eigen- 
functions shown in Figure[5l The constant c n is the ratio (g(x), 1}^ /\\g(x)\\fj 
for the n th eigenfunction. 



The resulting predictions for the Euler numbers are shown in the last column of Table [TJ 



5 Concluding remarks 



Another graph to investigate is the product with the even cycle C2 m , that is, £(Go s C2 m )- We 
conjecture that the resulting Euler number is asymptotically a constant times the associated Euler 
number for the product with a path, that is, 

8{Gn s C 2m ) _ 
£{Gu s P 2m ) ~^ c ' 
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Figure 6: A bipartite graph H m obtained by gluing m hexagons. 



as m tends to infinity and c is a positive constant less than 1 for S non-empty. 

Does the eigenfunction ip corresponding to the largest eigenvalue A carry information about the 
distribution of entries in the first copy of G in an alternating labeling of GO s P m ? More specifically, in 
the case of alternating 2 by m arrays, let £ (P2^Pm\ h j) be the number of alternating arrays where the 
first column consists of the two entries i and j, where 1 < i < j < 2m. Is the integer f (i^DP m ; i, j) 
approximated by c • (2m)! • A m_1 • (g(i/2m) + g(l — j /2m)) where c is the appropriate constant and 
g is the first eigenfunction displayed in Figure [5]? 

These techniques for obtaining the asymptotic behavior of the Euler numbers can be used for other 
classes of graphs as well. See for instance the graph H m in Figured! which is built by gluing hexagons 
together. Although Theorem 11.11 does not directly apply to this class of graphs, one can extend the 
theory to obtain the same asymptotic result. Hence we have 

£{Hm) ,m_l , nt ..m— 1\ 



c • A m + 0(//' 



(4-m + 2)! 

The essential question remaining is can the associated eigenvalue problem be solved explicitly. 

Keeping n fixed we know that £(P n nP m ) ~ C( n ) • (n • m)\ ■ A^" 1 for a constant C( n ) and the largest 
eigenvalue A( n ). Can anything be determined about the sequence A( n )? What can be said about the 
asymptotics of the Euler number £(P m OP m ) as m tends to infinity? 

A different direction is to study the descent number of directed graphs (digraphs). For a digraph 
G = (V, E) on n vertices define its descent number to be the number of labelings ir of the vertices with 
1 through n such that for each directed edge M-mwe have that tt(u) < n(v). If the digraph contains 
a directed cycle then the descent number is zero. For an acyclic digraph (digraphs without directed 
cycles) the descent number is strictly positive. The classical descent set statistics for permutations 
is obtained be looking at orientations of the path. By gluing directed graphs together, one obtains 
classes of graphs whose asymptotics of the descent number is natural to study via linear operators 
and their eigenvalues. 

The technique of translating a combinatorial problem into a problem of studying an operator and 
its spectrum was also explored in [1], where consecutive pattern avoiding in permutations were studied. 

Finally, we end with a purely enumerative question for trees (connected graphs without cycles). 
Conjecture 5.1. For a tree T on n vertices the classical Euler number E n is a lower bound for £ (T) , 
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that is, 

£{T) > E n . 

Furthermore, equality only holds when the tree T is the path P n . 
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